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ABSTRACT 



Context. Recent calculations have shown that the UV bump at about 217.5 nm in the extinction curve can be explained by a com- 
plex mixture of polycyclic aromatic hydrocarbons (PAHs) in several ionisation states. Other studies proposed that the carriers are a 
restricted population made of neutral and singly-ionised dehydrogenated coronene molecules (C24H„, n < 3), in line with models of 
the hydrogenation state of interstellar PAHs predicting that medium-sized species are highly dehydrogenated. 
Aims. To assess the observational consequences of the latter hypothesis we have undertaken a systematic theoretical study of the 
electronic spectra of dehydrogenated PAHs. We use our first results to see whether such spectra show strong general trends upon 
dehydrogenation. 

Methods. We performed calculations using state-of-the-art techniques in the framework of the density functional theory (DFT) to 
obtain the electronic ground-state geometries, and of the time-dependent DFT to evaluate the electronic excited-state properties. 
Results. We computed the absorption cross-section of the species C24H;^ (n= 12, 10,8,6,4,2,0) in their neutral and cationic 
charge-states. Similar calculations were performed for other PAHs and their fully dehydrogenated counterparts. 
Conclusions. ;7r-electron energies are always found to be strongly aff'ected by dehydrogenation. In all cases we examined, progressive 
dehydrogenation translates into a correspondingly progressive blue shift of the main electronic transitions. In particular, the tt ^ tt* 
collective resonance becomes broader and bluer with dehydrogenation. Its calculated energy position is therefore predicted to fall in 
the gap between the UV bump and the far-UV rise of the extinction curve. Since this eff'ect appears to be systematic, it poses a tight 
observational limit on the column density of strongly dehydrogenated medium-sized PAHs. 

Key words. (ISM:) dust, extinction -ISM: lines and bands -ISM: molecules -ultraviolet: ISM -astrochemistry -molecular data 



1. Introduction 



First discovered by lStecheil (Il965b , the absorption band centered 
at about 217.5 nm (~5.7 eV) is the most intense interstellar ex- 
tinction feature. With a relatively constant peak wavelen gth, the 
so-called "UV bump" has a Lorentz-like profile (e.g.. ISeatonl 
197 9) whose width varies considerably from one line of sight to 
another ( Eitzpatrick & Massa 2007). 

Both the strength and the position of the 217.5 nm fea- 
ture suggest that it originates in some form of carbonaceous 
material, known to display strong n ^ electron transi- 
tions in this range. A large number of carrier candidates have 
been proposed over the years, amongst which size-restricted 
graphite part icles ([Ste cher & Donn 1965 |), fuUerenes , C6o in par- 
ticular, (e.g., Braga et al. 1991; Iglesias-Groth 2004), polycyclic 
aromatic hydrocarbons (PAHs , Job lin 1992; Joblin et al. 1992), 
coal-like material (iPapoular et al.l Il996|). UV-processed hy - 
drogenated amorphous carbon particles ("Mennella et al.""l998), 
nano-sized hydroge nated carbon gr ains ( Schnaiter et al. 1998), 
and carbon onions ( C hhowalla et al.ll2QQ3l) . 

Since each PAH displays strong n transition in 

the 5 . 7 eV region 



119921: iJoblin et al.l Il992l : iDraind 

120031: iMalloci et al.1 [20041 ), mixtur es of PAH s must contribute 
to the 217.5 nm bump. Recently, C ecchi-P estellini et al. (2008) 
showed that such mixtures can account for the entire bump, as 
well as its strength relative to the non-linear far-UV rise, de- 



pending on the distribution of charge states in the PAH mix- 
tures. The UV bump was also attributed to a tt plas- 
mon resonance from a very limited number of strongly dehy- 
drogenated derivatives of coronene (C24H12), namely C24H„ and 
C24H„^ with n < 3. This assignment was made on the basis of 
DDA scattering calculations (Duley & Seahra 1998), supported 
by laboratory studies of the UV absorption in amorphous carbon 
(iDulev & LazarevI [2004) . The same molecules were als o pro- 
posed as carriers of some of the diflTuse interstellar bands (iDuleyl 
"2006b). Indeed, models of the hydrogenation state of PAHs pre- 
dict small and medium-siz ed PAHs (-20-30 carbon atoms) to 
be highly dehydrogenated ( Le Page et al.ll2QQ3b . 

Since the efifect of dehydrogenation on the electronic spectra 
of isolated PAHs is poorly characterized, we recently undertook 
a detailed study of the electronic excitation properties of de- 
hydro genated PAHs, using state-of-the-art quantum-chemical 
tools (iMalloci et al.l l2QQ8l) . As a part of this more detailed 
work, we present in this paper our first results for the dehydro- 
genated derivatives of coronene (C24H12), and completely de- 
hydrogenated perylene (C20H10) and ovalene (C32H14). We use 
them to asses s the specific case of the species introduced by 
[Puleyl (l2QQ6al) . whose computed spectral properties appear to be 
incompatible with the observed UV-bump, and search for gen- 
eral trends in the electron properties of dehydrogenated PAHs. 

Technical details of the calculations, as well as the justifica- 
tion for the choice of molecules considered, are given in Sect. [2 



2 
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Results are presented and discussed in the astrophysical context 
in Sect. O and our conclusions are reported in Sect. (H 



2. Methods 



2.1. Computational details 



As in previous studies (iMalloci et al.ll2QQ7allbb we used density 
functional theory (DFT Oones & Gunnarssonlll98 9) for the elec- 
tron ground state, and it s time-dependent extension (TD-DFT, 
iMarques & GrossI [2QQ4I) for electron excited states. Geometry 
optimisations were first performed using the computationally 
inexpensive 4-3 IG basis set, then refined with the 6-3 1+G* 
basis jFrish et al.lll984l) . F or this part o f the work we used the 
hybrid B3LYP functional (lBeckdll993h . as implement ed in the 
Gaussian-based DFT module of NWChem (lApra et al.ll20Q5l) . 

Starting from the geometry of each fully hydrogenated PAH 
(iMalloci et al. 2007a), we considered its dehydrogenated deriva- 
tives. We considered only neutral and cationic charge-states, 
and found in all cases that the computed ground-states for neu- 
trals and cations are the singlet and the doublet, respectively. 
All calculations were performed without symmetry constraints. 
To evaluate the eff'ect of dehydrogenation on the electronic 
structure of the species considered, we performed exploratory 
Harth ree-Fock/4-3 1 G calculations. Following Koopman's the- 
orem (lElandlll9^4l) . we interpreted the one-electron molecular 
orbital energies thus obtained as vertical ionisation energies. 

Keeping fixed the ground-state optimised geometries ob- 
tained above, we then computed the absorption cross-section 
cr(E) using the real-time TD-DFT implementation of octo- 
pus (Marques et al. 2003), in the local-density approximation 
dPerdew & Zungeii l 19811) as described in Malloci et al. (2007a). 
We previously showed (IMalloci et al.) 12004: Mulas et al. 2006) 
that the results of this method reproduce the overall far-UV be- 
haviour of cr(E), as seen in comparison with the experimental 
data available for a few neutral PAHs ( Joblin 1992; Joblin et al. 
[l992h . This includes the broad absorption peak dominated by 




(T ^ (T* transitions, which matches well both in position and 
width. The n ^ transitions from the visible to the UV range, 
while their absolute intensities are well reproduced, show ener- 
gies that are systematically underestimated. 

The TD-DFT formalism used in this work yields directly 
the time-dependent linear response of a given molecule after 
an impulsive perturbation, producing the whole spectrum in a 
single step (Marques et al. 2003). This diff'ers from the most 
widely used frequency-space TD-DFT implementation, where 
the poles of the linear response function correspond to vertical 
excitation energies, and the pol e strengths to the corresponding 
oscillator strengths (Casida 1995). We compare in Fig. [2] th e ex- 
perimental absorption s pectra of neu t ral naphthalene (Gingelll 



Il998h . and anthracene (iJoblinI Il992l: I Joblin 



. ^ )linetal.]ll992ir with 

both the real-time TD-DFT fMalloci et al.""2004') and the fre- 
quency-space TD-DFT (Malloci et al. 2007b) spectra. The lat- 
ter results agree better with the laboratory data. However, com- 
putational costs scale steeply with the number of transitions and 
electronic excitations of molecules as large as those considered 
here are thus limited to the visible/near-UV part of the spectrum. 




5 6 7 

Energy (eV) 



Fig. 2. Comparison between the experimental (black line) ab- 
sorption cross-section o'(E ) in me^abarn s (1Mb = 10"^^ cm^) 
of neutral naphthalene (top, Gingell 1998), and anthracene (bot- 
tom, Joblin 1992; Joblin et al. 1992), and two diff'erent the- 
oretical results: th e real-time TD-DFT spectrum (red line, 
IMalloci et al.l2004j). and t he frequency space TD-DFT one (blue 
line. IMalloci et al.ll2007bh . 



Fig. 1. Geometries of the dehydrogenated coronene molecules 
considered (C24H10, C24H8, C24H6, C24H4, C24H2, and C24), and 
of the completely dehydrogenated perylene (C20) and ovalene 
(C32). 



2.2. Choice of molecules 

We studied the dehydrogenated derivatives of coronene C24H„ 
(n= 10,8,6,4,2,0) in their neutral and cationic charge states. We 
included only species with an even n umber of hydrog en atoms, 
based on experimental mass spectra (Ekem et al.ll 19971) showing 



Malloci et al: Dehydrogenated PAHs and UV bump (RN) 



3 



that dehydrogenation corresponding to loss of an even number 
of hydrogen atoms is dominant. In addition, experimen tal re- 
sults on the dehydrogenation of the coronene cation (Joblin 'etaD 
bOQSh show that, for C24H2^ molecules, only species containing 
adjacent paired hy drogens have to be considered. In particular, 
iJoblin et aP (l2008b found that the main hydrogen loss channel is 
the detachment of single atoms; however, the photodissociation 
rate for "lone" hydrogen atoms was observed to be much faster 
(about one order of magnitude) than for paired hydrogen atoms. 
This implies that dehydrogenated coronene derivatives contain- 
ing adjacent paired hydrogens are more likely to be found in 
interstellar conditions. Moreover, it appears that the placement 
of the hydrogens at the periphery of the carbon skeleton has a 
minor effect on the main n ^ transitions. Such transitions 
depend mainly on the number of electrons populating the reso- 
nant TT and TT* molecular orbitals, which depend only on the num- 
ber of peripheral hydrogen atoms and not on their specific posi- 
tion. As an example. Fig. [3] shows the comparison between the 
computed absorption cross-sections of two isomers of C24Hio^, 
whose geometries are sketched in the same figure. As expected, 
with the exception of small details, the spectra of the two iso- 
mers are found to be almost indistinguishable in the energy range 
of astrophysical interest. Our study can therefore be restricted, 
without loss of generality for the scope of the present paper, to 
only the isomers of C24H„ (n= 10,8,6,4,2,0) shown in Fig.[T] We 
additionally examined the completely dehydrogenated deriva- 
tives of perylene (C20H10) and ovalene (C32H14) cations, i. e. 

and C32, representative of non-symmetric, dehydrogenated, 
compact, medium-sized PAHs. 



10 r 
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Fig. 3. Comparison between the computed absorption 
cross-section of the two inequivalent isomers of C24Hio^ 
sketched on the right. 



3. Results and discussion 

3.1. Dehydrogenation and electronic excitation properties 

According to Koopman's theorem the molecular orbital ener- 
gies obtained at the Harthree-Fock / 4-3 IG level are interpreted 
as vertical ionisation energies. We report in Fig. |4] the energies 
corresponding to removal of an electron from both the high- 
est occupied molecular orbital (IPl) and the molecular orbital 
just below it (IP2), as a function of the degree of hydrogena- 
tion (x = 0, 2, 4, 6, 8, 10, 12). In the model put forward by DulevI 
(l2006ah the low- and high-frequency cutoff's extracted from the 
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Fig. 4. First (asterisks) and second (diamonds) vertical ionisation 
energies of dehydrogenated coronene molecules C24B^x, com- 
puted at the Harthree-Fock/4-31G level in the framework of 
Koopman's theorem, as a function of the degree of hydrogena- 
tion. The horizontal dotted lines correspond to the low- and 
high-frequency cutoff's extracted fro m the interstellar extinction 
curve in the model bv iDuleyl (l2006ah . 



interstellar extinction curve, at about 7.3 and 9. 1 eV respectively, 
are identified with the IPl and IP2 of the carrier of the UV 
bump. These values were found to be in close agreement with 
those measured for ne utral coronene, 7.29 and 9.13 eV, respec- 
tively (IClar et al Based on the assumption that 7r-electron 
energies are little aff'ected by dehydrogenation, this supported 
the proposal that sl n ^ n* plasmon resonance in neutral and 
singly-ionised dehydrogenated coronene molecules could gen- 
erate the UV bump. However, we found an approximately con- 
stant change in ionization energy per hydrogen atom of about 
0.1 eV (cf. Fig. lU), suggesting that 7r-electron energies are in- 
deed strongly aff'ected by dehydrogenation. This is confirmed for 
the adiabatic first ionization energies computed via total energy 
diff'erences at the B3LYP/6-31-hG* level: 7.02, 7.23, 7.38, 7.61, 
7.79, 8.00, and 8.25 eV, when going from C24H12 to C24. 

The comparison between the computed absorption 
cross-sections of the diff'erent dehydrogenated coronene 
molecules considered, in their neutral and cationic 
charge-states, is presented in Fig. [S] The corresponding 
comparisons between C2oH[q and C32H[4 and their fully 
dehydrogenated counterparts are shown in Fig. [6j For both 
neutral and singly-ionised coronene derivatives, progressive 
dehydrogenation translates into a correspondingly progressive 
blue shift of the main electronic transitions (Fig. [5]). The n ^ 
collective resonance, in particular, becomes broader as well 
as bluer with dehydrogenation. The same holds true for the 
other PAHs reported in this work, namely perylene and ovalene 
(Fig. O, suggesting that this is a general feature of the whole 
class of molecules. 



3.2. Astrophysical implications 

For a quantitative comparison with extinction, we estimate the 
expected spectral contrast of the n ^ collective resonance 
in a mixture of all the fully dehydrogenated species considered 
so far. The top panel of Fig. [7] shows the weighted sum of their 
absorption spectra, compared with the same weighted sum for 
their fully hydrogenated parent molecules. The n ^ feature 
appears to peak at -6.95 eV (-178 nm), with an integrated in- 
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Fig. 5. Comparison between the computed absorption 
cross-sections of coronene C24H12 (black), C24H10 (violet), 
C24H8 (blue), C24H6 (cyan), C24H4 (green), C24H2 (orange), and 
C24 (red). Neutral molecules and their corresponding cations are 
displayed in the top and bottom panel, respectively. The dotted 
vertical line represents the position of the observed UV bump at 
-5.7 eV. 



tensity above the underlying continuum of about 1.3 times larger 
than the strong n ^ feature at 5.95 eV associated with the 
hydrogenated counterparts. 

Cecchi-Pestellini et al.' (l2008l) showed that mixtures of a 
large number of fully hydrogenated PAHs can reproduce the 
UV bump and non-linear far-UV rise of the extinction curve. 
In particular, the UV bump is very accurately matched in posi- 
tion (~5.7 eV), intensity and width. Such a fit can only be ob- 
tained using mixtures including a much larger sample of PAHs 
than those included here. However, the position of the n ^ 
resonance, and its intensity per carbon atom, are relatively insen- 
sitive to the individual species considered. In all cases, its peak is 
well within the observed width of the UV bump. This is the rea- 
son why such transitions, in a mixture of hydrogenated PAHs, 
become blended, yielding a smooth feature which can be very 
similar to the observed UV bump profile. 

In contrast, all completely dehydrogenated species consid- 
ered here have the peak of their n ^ resonance about ~ 1 eV 
bluewards of the bump, in the gap between it and the onset of the 
non-linear far-UV rise. The position and intensity per carbon 
atom of the n ^ resonances of dehydrogenated PAHs, as a 
family, can be expected to be as insensitive to individual species 
as for their fully hydrogenated counterparts. A larger mixture, if 
this is the case, would yield a smoother feature than the one we 
calculated with the present sample, but generally with the same 
intensity and position. 

Since no structure at ^ 7 eV is detected in the average galac- 
tic inte rstellar extinction curve reported by lFitzpatrick & Massal 
(l2007h . fully dehydrogenated PAHs must produce a peak smaller 



Fig. 6. Comparison between the computed absorption 
cross-sections of fully hydrogenated (black), and fully de- 
hydrogenated (red) perylene cation (C2oH[q) (top panel) and 
ovalene cation (C32H[4) (bottom panel). 



than the error of the interstellar extinction curve as illustrated 
in Fig. [71 Using a 2cr upper limit, we can derive that the inte- 
grated intensity of the n ^ feature of dehydrogenated PAHs 
is less than 1/6 of the bump area. Ass uming that the bump 
is m ainly due to hydrogenated PAHs ( Cecchi-Pestellini et al.l 
2008) and using the factor given above for the relative integrated 
cr{E), we derive that less than 1/8 of the total abundance of car- 
bon in PAHs is contained in strongly de hydro genated PAHs. 
Assuming Nh - 5.9 x 10^^ Eb-v (Whitte7"2002'), and applying 
the procedure of ICecchi^ estellini et al. (2008) t o the average 
extinction curve (Ry - 3.01 from Fitzpatrick & Massal l2007h 
yields a value of -150 ppM for the abundance of C locked in 
PAHs (including both the free-flying ones and those clustered 
in very small grains); the upper limit of the contribution of de- 
hydrogenated PAHs to the extinction curve can then be trans- 
lated into an upper limit of ~ 1 8 ppM for the carbon atoms con- 
tained in them. This numeric result is obviously dependent on 
the choice to compare against the average interstellar extinction 
curve, and would be difl'erent for a difl'erent choice. We showed 
in Cecchi-Pestellini et al. ( 2008) that the fraction of carbon con- 
tained in free-flying and clustered PAHs is quite variable, from 
85 to 187 ppm in a sample of lines of sight. 

Our up per limit can be compared with those found by 
IClayton et al.i (2003) regarding the column density of specific 
small, neutral, fully hydrogenated PAHs in specific lines of sight, 
based on laboratory spectra and Hubble Space Telescope obser- 
vations. Their results are more stringent thanks to the stronger 
spectral contrast of the absorption bands in their experimen- 
tal spectra, coupled with the matching high resolving power of 
the STIS instrument. For the present comparison, however, such 
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Fig. 7. Top panel: Weighted sum of the computed absorption 
cross-sections of C20, C24, C24, and C32 (continuous Hne), and 
C2oH[q, C24H12, C24H[2, and C32H[4 (dotted Hne), using as 
weight the inverse of the number of C atoms in each molecule. 
Bottom pane l: comparison between the average extinction curve 
(dotted line. iFitzpatrick & Massall2QQ7l ), with the Icr region in 
gray shade, and the average extinction curve + the amount of 
dehydrogenated PAHs which would make the n ^ feature 
detectable at 2cr (continuous line). 

higher resolving power would not help, given the broader nature 
of the spectral signature being sought. 

4. Conclusions 

PAHs are expected to exist in a wide variety of interstellar envi- 
ronments, in a complex statistical equ ilibriurn of diff erent charge 
and hydrogenation states (see e.g.. iTielensI l2Q05h . Modelling 
studies suggest that intermediate-size PAHs in the range of 
20-30 carbon atoms are stripped of most of their hydrogen atoms 
but still survive under the co nditions of the diffuse interstellar 
medium (iLe Page et al.ll2003h . 

We undertook a systematic theoretical study of the effects 
of dehydrogenation on the electron properties of neutral and 
ionised PAHs using state-of-the-art quantum-chemical tech- 
niques in the framework of the real-time TD-DFT. For all 
of the species considered so far, we found that, regardless of 
charge-state and molecular size, the n ^ collective reso- 
nance broadens and shifts to higher energies with increasing de- 
hydrogenation. 

When the molecule is half dehydrogenated or more, the band 
is shifted enough to fall between the bump and the far-UV rise 
of the extinction curve. The associated specie s cannot be the car- 
riers of the bump, contrary to the proposal of lDuleyl (l2006a). 

While it is overreaching to draw definitive conclusions based 
on the few cases considered, our upper limits are due to sys- 
tematic effects, and are thus unlikely to depend strongly on the 



specific sample of molec ules considered. No netheless, a larger 
study, which is underway ( Malloci et al. 2008), will be needed to 
assess to what extent the effect of dehydrogenation that we ob- 
served in UV spectra of PAHs is systematic for the whole class. 
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